knitr::opts_chunk$set(echo = TRUE)
rm(list = ls())
#cmdstanr::set_cmdstan_path(path = "C:/Users/kueng/.cmdstan/cmdstan-2.35.0")
#cmdstanr::set_cmdstan_path(path = "C:/Users/pascku/.cmdstan/cmdstan-2.36.0")
library(tidyverse)
library(R.utils)
library(wbCorr)
library(readxl)
library(kableExtra)
library(brms)
library(bayesplot)
library(see)
library(beepr)
library(DHARMa)
library(digest)
source(file.path('Functions', 'ReportModels.R'))
source(file.path('Functions', 'PrettyTables.R'))
source(file.path('Functions', 'ReportMeasures.R'))
source(file.path('Functions', 'PrepareData.R'))
report_function_hash <- digest::digest(summarize_brms)## [1] 1116
# Set options for analysis
use_mi = FALSE
shutdown = FALSE
report_ordinal = FALSE
do_priorsense = FALSE
get_bayesfactor = TRUE
check_models = TRUE #
if (get_bayesfactor) {
stats_to_report <- c('CI', 'SE', 'pd', 'ROPE', 'BF', 'Rhat', 'ESS')
} else {
stats_to_report <- c('CI', 'SE', 'pd', 'ROPE', 'Rhat', 'ESS')
}
options(
dplyr.print_max = 100,
brms.backend = 'cmdstan',
brms.file_refit = ifelse(use_mi, 'never', 'on_change'),
brms.file_refit = 'on_change',
#brms.file_refit = 'always',
error = function() {
beepr::beep(sound = 5)
if (shutdown) {
system("shutdown /s /t 180")
quit(save = "no", status = 1)
}
}
, es.use_symbols = TRUE
)
####################### Model parameters #######################
iterations = 12000 # 12'000 per chain to achieve 40'000
warmup = 2000 # 2000
# NO AR!!!
#corstr = 'ar'
#corstr = 'cosy_couple'
#corstr = 'cosy_couple:user'
################################################################
suffix = paste0('_sensitivityOtherCovariates_', as.character(iterations))df <- openxlsx::read.xlsx(file.path('long.xlsx'))
df_original <- df
df_double <- prepare_data(df, recode_pushing = TRUE, use_mi = use_mi)[[1]]Constructing scales Re-coding pusing reshaping data (4field) centering data within and between
Min. 1st Qu. Median Mean 3rd Qu. Max. NA’s 0.0000 0.0000 0.0000 0.1649 0.0000 5.0000 275
# For indistinguishable Dyads
model_rows_fixed <- c(
'Intercept',
# '-- WITHIN PERSON MAIN EFFECTS --',
'persuasion_self_cw',
'persuasion_partner_cw',
'pressure_self_cw',
'pressure_partner_cw',
'pushing_self_cw',
'pushing_partner_cw',
'day',
'plan_selfPlan',
'plan_partnerPlan',
'weartime_self_cw',
'support_self_cw',
'support_partner_cw',
'isWeekendWeekend',
'got_JITAI_selfJITAIreceived',
'skilled_supportDaysafterIntervention',
# '-- BETWEEN PERSON MAIN EFFECTS',
'persuasion_self_cb',
'persuasion_partner_cb',
'pressure_self_cb',
'pressure_partner_cb',
'pushing_self_cb',
'pushing_partner_cb',
'weartime_self_cb',
'support_self_cb',
'support_partner_cb',
'studyGroupFirst3weeksinterventions',
'studyGrouplast3weeksinterventions'
)
model_rows_fixed_ordinal <- c(
model_rows_fixed[1],
'Intercept[1]',
'Intercept[2]',
'Intercept[3]',
'Intercept[4]',
'Intercept[5]',
model_rows_fixed[2:length(model_rows_fixed)]
)
model_rows_random <- c(
# '--------------',
# '-- RANDOM EFFECTS --',
'sd(Intercept)',
'sd(persuasion_self_cw)',
'sd(persuasion_partner_cw)',
'sd(pressure_self_cw)',
'sd(pressure_partner_cw)',
'sd(pushing_self_cw)',
'sd(pushing_partner_cw)',
# '-- CORRELATION STRUCTURE -- ',
'sigma'
)
model_rows_random_ordinal <- c(model_rows_random,'disc')# For indistinguishable Dyads
model_rownames_fixed <- c(
"Intercept",
# "-- WITHIN PERSON MAIN EFFECTS --",
"Daily individual's experienced persuasion",
"Daily partner's experienced persuasion",
"Daily individual's experienced pressure",
"Daily partner's experienced pressure",
"Daily individual's experienced pushing",
"Daily partner's experienced pushing",
"Day",
"Own Actionplan",
'Partner Actionplan',
"Daily wear time",
'Daily support received',
"Daily support provided (partner's view)",
'Is a Weekend',
'JITAI received',
'Days post skilled support intervention',
# "-- BETWEEN PERSON MAIN EFFECTS",
"Mean individual's experienced persuasion",
"Mean partner's experienced persuasion",
"Mean individual's experienced pressure",
"Mean partner's experienced pressure",
"Mean individual's experienced pushing",
"Mean partner's experienced pushing",
"Mean wear time",
"Mean individual's received support",
"Mean partner's received support",
'Difference study group 2',
'Difference study group 3'
)
model_rownames_fixed_ordinal <- c(
model_rownames_fixed[1],
'Intercept[1]',
'Intercept[2]',
'Intercept[3]',
'Intercept[4]',
'Intercept[5]',
model_rownames_fixed[2:length(model_rownames_fixed)]
)
model_rownames_random <- c(
# '--------------',
# '-- RANDOM EFFECTS --',
'sd(Intercept)',
"sd(Daily individual's experienced persuasion)",
"sd(Daily partner's experienced persuasion)", # OR partner received
"sd(Daily individual's experienced pressure)",
"sd(Daily partner's experienced pressure)",
"sd(Daily individual's experienced pushing)",
"sd(Daily partner's experienced pushing)",
# '-- CORRELATION STRUCTURE -- ',
'sigma'
)
model_rownames_random_ordinal <- c(model_rownames_random,'disc')rows_to_pack <- list(
"Within-Person Effects" = c(2,16),
"Between-Person Effects" = c(17,27),
"Random Effects" = c(28, 34),
"Additional Parameters" = c(35,35)
)
rows_to_pack_ordinal <- list(
"Intercepts" = c(1,6),
"Within-Person Effects" = c(2+5,16+5),
"Between-Person Effects" = c(17+5,27+5),
"Random Effects" = c(28+5, 34+5),
"Additional Parameters" = c(35+5,35+6)
)formula <- bf(
pa_sub ~
persuasion_self_cw + persuasion_partner_cw +
pressure_self_cw + pressure_partner_cw +
pushing_self_cw + pushing_partner_cw +
#support_self_cw + support_partner_cw +
isWeekend +
got_JITAI_self +
skilled_support +
persuasion_self_cb + persuasion_partner_cb +
pressure_self_cb + pressure_partner_cb +
pushing_self_cb + pushing_partner_cb +
#support_self_cb + support_partner_cb +
studyGroup +
plan_self + plan_partner +
day +
# Random effects
(1 + persuasion_self_cw + persuasion_partner_cw +
pressure_self_cw + pressure_partner_cw +
pushing_self_cw + pushing_partner_cw | dd | coupleID),
hu = ~ persuasion_self_cw + persuasion_partner_cw +
pressure_self_cw + pressure_partner_cw +
pushing_self_cw + pushing_partner_cw +
#support_self_cw + support_partner_cw +
isWeekend +
got_JITAI_self +
skilled_support +
persuasion_self_cb + persuasion_partner_cb +
pressure_self_cb + pressure_partner_cb +
pushing_self_cb + pushing_partner_cb +
#support_self_cb + support_partner_cb +
studyGroup +
plan_self + plan_partner +
day +
# Random effects
(1 + persuasion_self_cw + persuasion_partner_cw +
pressure_self_cw + pressure_partner_cw +
pushing_self_cw + pushing_partner_cw | dd | coupleID)
, decomp = 'QR'
#, autocor = autocor_str
)
prior1 <- c(
brms::set_prior("normal(0, 2.5)", class = "b")
, brms::set_prior("normal(0, 2)", class = "b", dpar = "hu")
, brms::set_prior("normal(0, 50)", class = "Intercept") # for non-zero PA
, brms::set_prior("normal(0.5, 2.5)", class = "Intercept", dpar = 'hu') # hurdle part
, brms::set_prior("normal(0, 2)", class = "sd", group = "coupleID", lb = 0)
, brms::set_prior("student_t(3, 0, 2.5)", class = "sigma", lb = 0)
)
#brms::validate_prior(
# prior1,
# formula = formula,
# data = df_double,
# family = hurdle_lognormal()
#)
#df_minimal <- df_double[, c("AorB", all.vars(as.formula(formula)))]
pa_sub <- my_brm(
mi = use_mi,
imputed_data = implist,
formula = formula,
prior = prior1,
data = df_double,
family = brms::hurdle_lognormal(),
#family = brms::hurdle_negbinomial(),
#family = brms::hurdle_poisson(),
#control = list(adapt_delta = 0.95, max_treedepth = 15),
iter = iterations,
warmup = warmup,
chains = 4,
cores = 4,
seed = 42,
file = file.path("models_cache_brms", paste0("pa_sub_hu_lognormal", suffix))
#, file_refit = 'always'
)## Warning: Rows containing NAs were excluded from the model.
# rope range for continuous part of the model
rope_factor <- sd(log(pa_sub$data$pa_sub[pa_sub$data$pa_sub > 0]))
rope_range_continuous = c(-0.1 * rope_factor, 0.1 * rope_factor)
summary_pa_sub <- summarize_brms(
pa_sub,
stats_to_report = stats_to_report,
rope_range = rope_range_continuous,
hu_rope_range = c(-0.18, 0.18),
model_rows_fixed = model_rows_fixed,
model_rows_random = model_rows_random,
model_rownames_fixed = model_rownames_fixed,
model_rownames_random = model_rownames_random,
exponentiate = T) ## Sampling priors, please wait...
## Warning in summarize_brms(pa_sub, stats_to_report = stats_to_report,
## rope_range = rope_range_continuous, : Coefficients were exponentiated.
## Double check if this was intended.
| exp(Est.)_hu | SE_hu | 95% CI_hu | pd_hu | ROPE_hu | inside ROPE_hu | BF_hu | BF_Evidence_hu | Rhat_hu | Bulk_ESS_hu | Tail_ESS_hu | exp(Est.)_nonzero | SE_nonzero | 95% CI_nonzero | pd_nonzero | ROPE_nonzero | inside ROPE_nonzero | BF_nonzero | BF_Evidence_nonzero | Rhat_nonzero | Bulk_ESS_nonzero | Tail_ESS_nonzero | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Intercept | 0.18*** | 0.05 | [ 0.10, 0.30] | 1.000 | [0.84, 1.20] | 0.000 | >100 | Overwhelming Evidence | 1.000 | 35733 | 31253 | 29.02*** | 3.33 | [23.19, 36.50] | 1.000 | [0.92, 1.08] | 0.000 | >100 | Overwhelming Evidence | 1.000 | 18775 | 26518 |
| Within-Person Effects | ||||||||||||||||||||||
| Daily individual’s experienced persuasion | 1.57*** | 0.11 | [ 1.38, 1.83] | 1.000 | [0.84, 1.20] | 0.000 | >100 | Overwhelming Evidence | 1.000 | 51560 | 33810 | 1.03 | 0.03 | [ 0.98, 1.08] | 0.855 | [0.92, 1.08] | 0.967 | 0.016 | Very Strong Evidence for Null | 1.000 | 34176 | 33120 |
| Daily partner’s experienced persuasion | 1.33*** | 0.09 | [ 1.18, 1.53] | 1.000 | [0.84, 1.20] | 0.046 | >100 | Overwhelming Evidence | 1.000 | 51987 | 30034 | 1.03 | 0.02 | [ 0.98, 1.07] | 0.871 | [0.92, 1.08] | 0.986 | 0.016 | Very Strong Evidence for Null | 1.000 | 39777 | 34751 |
| Daily individual’s experienced pressure | 0.98 | 0.16 | [ 0.70, 1.37] | 0.540 | [0.84, 1.20] | 0.724 | 0.080 | Strong Evidence for Null | 1.000 | 53291 | 32379 | 0.90* | 0.05 | [ 0.81, 1.00] | 0.978 | [0.92, 1.08] | 0.293 | 0.066 | Strong Evidence for Null | 1.000 | 41024 | 37149 |
| Daily partner’s experienced pressure | 1.56* | 0.34 | [ 1.04, 2.74] | 0.983 | [0.84, 1.20] | 0.101 | 1.054 | Weak Evidence | 1.000 | 40719 | 27620 | 0.95 | 0.04 | [ 0.86, 1.04] | 0.880 | [0.92, 1.08] | 0.715 | 0.012 | Very Strong Evidence for Null | 1.000 | 42700 | 38100 |
| Daily individual’s experienced pushing | 1.00 | 0.14 | [ 0.76, 1.35] | 0.501 | [0.84, 1.20] | 0.785 | 0.071 | Strong Evidence for Null | 1.000 | 46932 | 33301 | 0.99 | 0.03 | [ 0.92, 1.05] | 0.662 | [0.92, 1.08] | 0.972 | 0.008 | Very Strong Evidence for Null | 1.000 | 49562 | 36055 |
| Daily partner’s experienced pushing | 1.33* | 0.16 | [ 1.05, 1.72] | 0.992 | [0.84, 1.20] | 0.189 | 1.125 | Weak Evidence | 1.000 | 55699 | 32118 | 0.96 | 0.03 | [ 0.90, 1.02] | 0.903 | [0.92, 1.08] | 0.898 | 0.015 | Very Strong Evidence for Null | 1.000 | 49205 | 33720 |
| Day | 0.69 | 0.15 | [ 0.46, 1.04] | 0.961 | [0.84, 1.20] | 0.178 | 0.487 | Weak Evidence for Null | 1.000 | 67971 | 32562 | 0.94 | 0.08 | [ 0.79, 1.11] | 0.763 | [0.92, 1.08] | 0.525 | 0.009 | Very Strong Evidence for Null | 1.000 | 71902 | 32082 |
| Own Actionplan | 11.49*** | 1.23 | [ 9.35, 14.19] | 1.000 | [0.84, 1.20] | 0.000 | >100 | Overwhelming Evidence | 1.000 | 71441 | 30644 | 1.37*** | 0.07 | [ 1.24, 1.50] | 1.000 | [0.92, 1.08] | 0.000 | >100 | Overwhelming Evidence | 1.000 | 83623 | 32212 |
| Partner Actionplan | 1.24* | 0.13 | [ 1.01, 1.51] | 0.980 | [0.84, 1.20] | 0.374 | 0.407 | Weak Evidence for Null | 1.000 | 79370 | 32484 | 1.08 | 0.05 | [ 0.99, 1.18] | 0.965 | [0.92, 1.08] | 0.496 | 0.038 | Strong Evidence for Null | 1.000 | 79643 | 31132 |
| Daily wear time | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Daily support received | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Daily support provided (partner’s view) | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Is a Weekend | 1.57*** | 0.16 | [ 1.29, 1.91] | 1.000 | [0.84, 1.20] | 0.003 | >100 | Overwhelming Evidence | 1.000 | 91266 | 30219 | 1.19*** | 0.05 | [ 1.10, 1.28] | 1.000 | [0.92, 1.08] | 0.007 | 29.126 | Strong Evidence | 1.000 | 86831 | 28268 |
| JITAI received | 0.49*** | 0.06 | [ 0.38, 0.63] | 1.000 | [0.84, 1.20] | 0.000 | >100 | Overwhelming Evidence | 1.000 | 82907 | 32045 | 0.94 | 0.05 | [ 0.85, 1.04] | 0.874 | [0.92, 1.08] | 0.632 | 0.013 | Very Strong Evidence for Null | 1.000 | 66812 | 30880 |
| Days post skilled support intervention | 1.31 | 0.21 | [ 0.96, 1.80] | 0.954 | [0.84, 1.20] | 0.281 | 0.329 | Weak Evidence for Null | 1.000 | 65743 | 31177 | 1.07 | 0.07 | [ 0.94, 1.21] | 0.853 | [0.92, 1.08] | 0.573 | 0.012 | Very Strong Evidence for Null | 1.000 | 55612 | 36199 |
| Between-Person Effects | ||||||||||||||||||||||
| Mean individual’s experienced persuasion | 1.44 | 0.57 | [ 0.65, 3.13] | 0.815 | [0.84, 1.20] | 0.233 | 0.297 | Moderate Evidence for Null | 1.000 | 28924 | 30569 | 0.99 | 0.15 | [ 0.72, 1.34] | 0.536 | [0.92, 1.08] | 0.390 | 0.014 | Very Strong Evidence for Null | 1.000 | 30611 | 31893 |
| Mean partner’s experienced persuasion | 1.32 | 0.52 | [ 0.60, 2.88] | 0.760 | [0.84, 1.20] | 0.275 | 0.254 | Moderate Evidence for Null | 1.000 | 28495 | 28933 | 0.93 | 0.14 | [ 0.68, 1.26] | 0.685 | [0.92, 1.08] | 0.352 | 0.016 | Very Strong Evidence for Null | 1.000 | 30014 | 30924 |
| Mean individual’s experienced pressure | 0.16*** | 0.09 | [ 0.05, 0.47] | 1.000 | [0.84, 1.20] | 0.001 | 70.857 | Very Strong Evidence | 1.000 | 40609 | 34634 | 0.69 | 0.26 | [ 0.33, 1.45] | 0.840 | [0.92, 1.08] | 0.100 | 0.037 | Strong Evidence for Null | 1.000 | 12525 | 19285 |
| Mean partner’s experienced pressure | 0.27* | 0.15 | [ 0.08, 0.78] | 0.992 | [0.84, 1.20] | 0.016 | 5.019 | Moderate Evidence | 1.000 | 40177 | 33679 | 0.53 | 0.20 | [ 0.25, 1.12] | 0.953 | [0.92, 1.08] | 0.039 | 0.091 | Strong Evidence for Null | 1.000 | 12126 | 18855 |
| Mean individual’s experienced pushing | 1.39 | 0.85 | [ 0.42, 4.73] | 0.708 | [0.84, 1.20] | 0.204 | 0.350 | Weak Evidence for Null | 1.000 | 36971 | 33082 | 1.47 | 0.40 | [ 0.86, 2.52] | 0.922 | [0.92, 1.08] | 0.082 | 0.042 | Strong Evidence for Null | 1.000 | 17225 | 24575 |
| Mean partner’s experienced pushing | 2.33 | 1.41 | [ 0.70, 7.85] | 0.918 | [0.84, 1.20] | 0.090 | 0.799 | Weak Evidence for Null | 1.000 | 36654 | 32185 | 1.63 | 0.44 | [ 0.94, 2.83] | 0.960 | [0.92, 1.08] | 0.048 | 0.073 | Strong Evidence for Null | 1.000 | 18057 | 23622 |
| Mean wear time | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Mean individual’s received support | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Mean partner’s received support | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Difference study group 2 | 0.79 | 0.28 | [ 0.39, 1.60] | 0.752 | [0.84, 1.20] | 0.316 | 0.224 | Moderate Evidence for Null | 1.000 | 33038 | 30920 | 1.06 | 0.14 | [ 0.81, 1.38] | 0.672 | [0.92, 1.08] | 0.412 | 0.023 | Very Strong Evidence for Null | 1.000 | 19899 | 26220 |
| Difference study group 3 | 1.17 | 0.39 | [ 0.59, 2.33] | 0.675 | [0.84, 1.20] | 0.369 | 0.186 | Moderate Evidence for Null | 1.000 | 35675 | 29635 | 1.25 | 0.16 | [ 0.97, 1.60] | 0.962 | [0.92, 1.08] | 0.110 | 0.091 | Strong Evidence for Null | 1.000 | 21415 | 25449 |
| Random Effects | ||||||||||||||||||||||
| sd(Intercept) | 0.83 | 0.12 | [0.63, 1.12] | NA | NA | NA | NA | NA | 1.000 | 26004 | 29358 | 0.29 | 0.04 | [0.22, 0.40] | NA | NA | NA | NA | NA | 1.000 | 18362 | 26998 |
| sd(Daily individual’s experienced persuasion) | 0.20 | 0.09 | [0.03, 0.40] | NA | NA | NA | NA | NA | 1.000 | 17340 | 18974 | 0.12 | 0.02 | [0.08, 0.17] | NA | NA | NA | NA | NA | 1.000 | 27247 | 31529 |
| sd(Daily partner’s experienced persuasion) | 0.16 | 0.10 | [0.01, 0.36] | NA | NA | NA | NA | NA | 1.000 | 14203 | 20274 | 0.08 | 0.02 | [0.04, 0.13] | NA | NA | NA | NA | NA | 1.000 | 27278 | 27450 |
| sd(Daily individual’s experienced pressure) | 0.21 | 0.19 | [0.01, 0.78] | NA | NA | NA | NA | NA | 1.000 | 25017 | 26936 | 0.07 | 0.06 | [0.00, 0.23] | NA | NA | NA | NA | NA | 1.000 | 22361 | 24734 |
| sd(Daily partner’s experienced pressure) | 0.39 | 0.32 | [0.02, 1.25] | NA | NA | NA | NA | NA | 1.000 | 18619 | 25027 | 0.05 | 0.05 | [0.00, 0.18] | NA | NA | NA | NA | NA | 1.000 | 27800 | 25034 |
| sd(Daily individual’s experienced pushing) | 0.51 | 0.17 | [0.21, 0.91] | NA | NA | NA | NA | NA | 1.000 | 24757 | 22151 | 0.08 | 0.04 | [0.01, 0.16] | NA | NA | NA | NA | NA | 1.000 | 16136 | 17692 |
| sd(Daily partner’s experienced pushing) | 0.23 | 0.16 | [0.01, 0.59] | NA | NA | NA | NA | NA | 1.000 | 18975 | 22686 | 0.07 | 0.04 | [0.01, 0.14] | NA | NA | NA | NA | NA | 1.000 | 17731 | 17041 |
| Additional Parameters | ||||||||||||||||||||||
| sigma | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | 0.67 | 0.01 | [0.65, 0.70] | NA | NA | NA | NA | NA | 1.000 | 67759 | 30466 |
# Plot continuous part of model
variable <- c(
'(Intercept)',
'b_persuasion_self_cw',
'b_persuasion_partner_cw',
'b_pressure_self_cw',
'b_pressure_partner_cw',
'b_pushing_self_cw',
'b_pushing_partner_cw'
)
plot(
bayestestR::p_direction(pa_sub, parameter = variable),
priors = TRUE
) + theme_bw()plot(
bayestestR::rope(
pa_sub,
parameter = variable,
range = rope_range_continuous,
verbose = F,
ci = 1
)
) + theme_bw()# Hurdle part of the model
variable <- c(
'b_hu_persuasion_self_cw',
'b_hu_persuasion_partner_cw',
'b_hu_pressure_self_cw',
'b_hu_pressure_partner_cw',
'b_hu_pushing_self_cw',
'b_hu_pushing_partner_cw'
)
plot(
bayestestR::p_direction(pa_sub, parameter = variable),
priors = TRUE
) + theme_bw()# The rope range for the bernoulli part of the model is -0.18, 0.18
plot(
bayestestR::rope(pa_sub, parameter = variable, range = c(-0.18, 0.18), ci = 1),
verbose = FALSE
) + theme_bw()## Possible multicollinearity between b_hu_persuasion_partner_cb and
## b_hu_persuasion_self_cb (r = 0.78), b_persuasion_partner_cb and
## b_persuasion_self_cb (r = 0.76), b_pressure_partner_cb and
## b_pressure_self_cb (r = 0.82). This might lead to inappropriate
## results. See 'Details' in '?rope'.
formula <- bf(
pa_obj ~
persuasion_self_cw + persuasion_partner_cw +
pressure_self_cw + pressure_partner_cw +
pushing_self_cw + pushing_partner_cw +
#support_self_cw + support_partner_cw +
isWeekend +
got_JITAI_self +
skilled_support +
persuasion_self_cb + persuasion_partner_cb +
pressure_self_cb + pressure_partner_cb +
pushing_self_cb + pushing_partner_cb +
#support_self_cb + support_partner_cb +
studyGroup +
plan_self + plan_partner +
day + + weartime_self_cw + weartime_self_cb +
# Random effects
(persuasion_self_cw + persuasion_partner_cw +
pressure_self_cw + pressure_partner_cw +
pushing_self_cw + pushing_partner_cw | coupleID)
, decomp = 'QR'
#, autocor = autocor_str
)
prior1 <- c(
brms::set_prior("normal(0, 2.5)", class = "b")
, brms::set_prior("normal(0, 50)", class = "Intercept")
, brms::set_prior("normal(0, 2)", class = "sd", group = "coupleID", lb = 0)
, brms::set_prior("student_t(3, 0, 2.5)", class = "sigma", lb = 0)
)
#brms::validate_prior(
# prior1,
# formula = formula,
# data = df_double,
# family = lognormal()
# )
#df_minimal <- df_double[, c("AorB", all.vars(as.formula(formula)))]
pa_obj_log <- my_brm(
mi = use_mi,
imputed_data = implist,
formula = formula,
prior = prior1,
data = df_double,
family = lognormal(),
#control = list(adapt_delta = 0.95),
iter = iterations,
warmup = warmup,
chains = 4,
cores = 4,
seed = 7777,
file = file.path("models_cache_brms", paste0("pa_obj_log_gaussian", suffix))
)## Warning: Rows containing NAs were excluded from the model.
# rope range for lognormal model
rope_factor <- sd(log(pa_obj_log$data$pa_obj))
rope_range_log = c(-0.1 * rope_factor, 0.1 * rope_factor)
summary_pa_obj <- summarize_brms(
pa_obj_log,
stats_to_report = stats_to_report,
rope_range = rope_range_log,
model_rows_fixed = model_rows_fixed,
model_rows_random = model_rows_random,
model_rownames_fixed = model_rownames_fixed,
model_rownames_random = model_rownames_random,
exponentiate = T) ## Warning in compute_rope(range = rope_range): Collinearity detected. Some
## VIFs are > 10. This may invalidate ROPE inferences!
Low Correlation
Term VIF VIF 95% CI Increased SE Tolerance
persuasion_self_cw 1.11 [ 1.08, 1.15] 1.05 0.90
persuasion_partner_cw 1.14 [ 1.10, 1.18] 1.07 0.88 pressure_self_cw 1.07 [ 1.04, 1.12] 1.04 0.93 pressure_partner_cw 1.11 [ 1.08, 1.16] 1.06 0.90 pushing_self_cw 1.13 [ 1.09, 1.17] 1.06 0.89 pushing_partner_cw 1.18 [ 1.14, 1.23] 1.09 0.85 isWeekend 1.05 [ 1.02, 1.10] 1.02 0.96 got_JITAI_self 1.09 [ 1.06, 1.13] 1.04 0.92 skilled_support 2.01 [ 1.92, 2.10] 1.42 0.50 pressure_self_cb 4.28 [ 4.05, 4.52] 2.07 0.23 pressure_partner_cb 4.69 [ 4.44, 4.96] 2.17 0.21 studyGroup 1.17 [ 1.13, 1.21] 1.08 0.86 plan_self 1.30 [ 1.25, 1.35] 1.14 0.77 plan_partner 1.30 [ 1.25, 1.35] 1.14 0.77 day 1.94 [ 1.85, 2.03] 1.39 0.52 weartime_self_cw 1.05 [ 1.03, 1.10] 1.03 0.95 weartime_self_cb 1.25 [ 1.21, 1.30] 1.12 0.80 Tolerance 95% CI [0.87, 0.93] [0.85, 0.91] [0.90, 0.96] [0.86, 0.93] [0.85, 0.92] [0.82, 0.88] [0.91, 0.98] [0.88, 0.95] [0.48, 0.52] [0.22, 0.25] [0.20, 0.23] [0.82, 0.88] [0.74, 0.80] [0.74, 0.80] [0.49, 0.54] [0.91, 0.97] [0.77, 0.83]
Moderate Correlation
Term VIF VIF 95% CI Increased SE Tolerance
persuasion_self_cb 9.97 [ 9.40, 10.57] 3.16 0.10 pushing_self_cb 6.75 [ 6.38, 7.15] 2.60 0.15 pushing_partner_cb 6.58 [ 6.22, 6.97] 2.57 0.15 Tolerance 95% CI [0.09, 0.11] [0.14, 0.16] [0.14, 0.16]
High Correlation
Term VIF VIF 95% CI Increased SE Tolerance
persuasion_partner_cb 10.84 [10.22, 11.49] 3.29 0.09 Tolerance 95% CI [0.09, 0.10]
## Sampling priors, please wait...
## Warning in summarize_brms(pa_obj_log, stats_to_report = stats_to_report, :
## Coefficients were exponentiated. Double check if this was intended.
| exp(Est.) | SE | 95% CI | pd | ROPE | inside ROPE | BF | BF_Evidence | Rhat | Bulk_ESS | Tail_ESS | |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Intercept | 108.21*** | 9.80 | [90.06, 129.44] | 1.000 | [0.94, 1.07] | 0.000 | >100 | Overwhelming Evidence | 1.000 | 8093 | 13915 |
| Within-Person Effects | |||||||||||
| Daily individual’s experienced persuasion | 1.03 | 0.02 | [ 1.00, 1.06] | 0.951 | [0.94, 1.07] | 0.996 | 0.021 | Very Strong Evidence for Null | 1.000 | 33778 | 30856 |
| Daily partner’s experienced persuasion | 1.01 | 0.02 | [ 0.98, 1.05] | 0.794 | [0.94, 1.07] | 0.999 | 0.008 | Very Strong Evidence for Null | 1.000 | 34387 | 29786 |
| Daily individual’s experienced pressure | 0.95 | 0.03 | [ 0.88, 1.01] | 0.945 | [0.94, 1.07] | 0.611 | 0.018 | Very Strong Evidence for Null | 1.000 | 52258 | 32119 |
| Daily partner’s experienced pressure | 0.98 | 0.03 | [ 0.92, 1.05] | 0.705 | [0.94, 1.07] | 0.909 | 0.005 | Very Strong Evidence for Null | 1.000 | 57383 | 31698 |
| Daily individual’s experienced pushing | 1.02 | 0.03 | [ 0.96, 1.07] | 0.728 | [0.94, 1.07] | 0.971 | 0.007 | Very Strong Evidence for Null | 1.000 | 35362 | 29553 |
| Daily partner’s experienced pushing | 1.00 | 0.02 | [ 0.96, 1.04] | 0.512 | [0.94, 1.07] | 0.997 | 0.004 | Very Strong Evidence for Null | 1.000 | 56192 | 30841 |
| Day | 0.94 | 0.05 | [ 0.86, 1.03] | 0.897 | [0.94, 1.07] | 0.523 | 0.009 | Very Strong Evidence for Null | 1.000 | 79678 | 29938 |
| Own Actionplan | 1.08** | 0.03 | [ 1.02, 1.13] | 0.998 | [0.94, 1.07] | 0.369 | 0.303 | Weak Evidence for Null | 1.000 | 81270 | 27656 |
| Partner Actionplan | 1.05* | 0.03 | [ 1.00, 1.11] | 0.983 | [0.94, 1.07] | 0.681 | 0.039 | Strong Evidence for Null | 1.000 | 75784 | 29403 |
| Daily wear time | 1.00*** | 0.00 | [ 1.00, 1.00] | 1.000 | [0.94, 1.07] | 1.000 | 33.548 | Very Strong Evidence | 1.000 | 84568 | 28829 |
| Daily support received | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Daily support provided (partner’s view) | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Is a Weekend | 1.06** | 0.02 | [ 1.01, 1.11] | 0.995 | [0.94, 1.07] | 0.579 | 0.117 | Moderate Evidence for Null | 1.000 | 72368 | 28911 |
| JITAI received | 0.93* | 0.03 | [ 0.87, 0.99] | 0.993 | [0.94, 1.07] | 0.380 | 0.078 | Strong Evidence for Null | 1.000 | 76607 | 31051 |
| Days post skilled support intervention | 1.05 | 0.04 | [ 0.98, 1.13] | 0.915 | [0.94, 1.07] | 0.653 | 0.010 | Very Strong Evidence for Null | 1.000 | 55394 | 35489 |
| Between-Person Effects | |||||||||||
| Mean individual’s experienced persuasion | 1.10 | 0.16 | [ 0.82, 1.49] | 0.740 | [0.94, 1.07] | 0.279 | 0.016 | Very Strong Evidence for Null | 1.000 | 8034 | 14915 |
| Mean partner’s experienced persuasion | 0.98 | 0.15 | [ 0.73, 1.33] | 0.546 | [0.94, 1.07] | 0.336 | 0.014 | Very Strong Evidence for Null | 1.000 | 8049 | 14874 |
| Mean individual’s experienced pressure | 0.98 | 0.15 | [ 0.72, 1.32] | 0.565 | [0.94, 1.07] | 0.324 | 0.008 | Very Strong Evidence for Null | 1.000 | 10867 | 19807 |
| Mean partner’s experienced pressure | 0.98 | 0.14 | [ 0.73, 1.32] | 0.549 | [0.94, 1.07] | 0.337 | 0.009 | Very Strong Evidence for Null | 1.000 | 9460 | 17590 |
| Mean individual’s experienced pushing | 0.96 | 0.20 | [ 0.63, 1.47] | 0.577 | [0.94, 1.07] | 0.236 | 0.011 | Very Strong Evidence for Null | 1.000 | 10569 | 17814 |
| Mean partner’s experienced pushing | 1.20 | 0.25 | [ 0.79, 1.84] | 0.812 | [0.94, 1.07] | 0.163 | 0.017 | Very Strong Evidence for Null | 1.000 | 10506 | 17487 |
| Mean wear time | 1.00 | 0.00 | [ 1.00, 1.00] | 0.934 | [0.94, 1.07] | 1.000 | 0.021 | Very Strong Evidence for Null | 1.000 | 43730 | 31904 |
| Mean individual’s received support | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Mean partner’s received support | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Difference study group 2 | 0.94 | 0.11 | [ 0.73, 1.20] | 0.704 | [0.94, 1.07] | 0.349 | 0.023 | Very Strong Evidence for Null | 1.000 | 7388 | 13983 |
| Difference study group 3 | 1.03 | 0.13 | [ 0.80, 1.33] | 0.604 | [0.94, 1.07] | 0.388 | 0.017 | Very Strong Evidence for Null | 1.001 | 7540 | 13272 |
| Random Effects | |||||||||||
| sd(Intercept) | 0.30 | 0.04 | [0.24, 0.40] | NA | NA | NA | NA | NA | 1.000 | 10270 | 18425 |
| sd(Daily individual’s experienced persuasion) | 0.05 | 0.01 | [0.02, 0.08] | NA | NA | NA | NA | NA | 1.000 | 18370 | 15467 |
| sd(Daily partner’s experienced persuasion) | 0.05 | 0.02 | [0.02, 0.09] | NA | NA | NA | NA | NA | 1.000 | 16543 | 11878 |
| sd(Daily individual’s experienced pressure) | 0.04 | 0.03 | [0.00, 0.13] | NA | NA | NA | NA | NA | 1.000 | 21491 | 21233 |
| sd(Daily partner’s experienced pressure) | 0.03 | 0.03 | [0.00, 0.11] | NA | NA | NA | NA | NA | 1.000 | 24541 | 19139 |
| sd(Daily individual’s experienced pushing) | 0.08 | 0.04 | [0.01, 0.16] | NA | NA | NA | NA | NA | 1.000 | 10016 | 12350 |
| sd(Daily partner’s experienced pushing) | 0.03 | 0.02 | [0.00, 0.09] | NA | NA | NA | NA | NA | 1.000 | 16746 | 18525 |
| Additional Parameters | |||||||||||
| sigma | 0.57 | 0.01 | [0.56, 0.59] | NA | NA | NA | NA | NA | 1.000 | 64506 | 27834 |
plot(
bayestestR::p_direction(pa_obj_log),
priors = TRUE
) +
coord_cartesian(xlim = c(-3, 3)) +
theme_bw()## Warning in `==.default`(dens$Parameter, parameter): longer object length
## is not a multiple of shorter object length
## Warning in is.na(e1) | is.na(e2): longer object length is not a multiple
## of shorter object length
## Possible multicollinearity between b_persuasion_partner_cb and
## b_persuasion_self_cb (r = 0.9), b_pressure_self_cb and
## b_persuasion_self_cb (r = 0.73), b_pressure_partner_cb and
## b_persuasion_self_cb (r = 0.74), b_pressure_self_cb and
## b_persuasion_partner_cb (r = 0.71), b_pressure_partner_cb and
## b_persuasion_partner_cb (r = 0.77), b_pushing_partner_cb and
## b_pushing_self_cb (r = 0.84). This might lead to inappropriate
## results. See 'Details' in '?rope'.
formula <- bf(
aff ~
persuasion_self_cw + persuasion_partner_cw +
pressure_self_cw + pressure_partner_cw +
pushing_self_cw + pushing_partner_cw +
#support_self_cw + support_partner_cw +
isWeekend +
got_JITAI_self +
skilled_support +
persuasion_self_cb + persuasion_partner_cb +
pressure_self_cb + pressure_partner_cb +
pushing_self_cb + pushing_partner_cb +
#support_self_cb + support_partner_cb +
studyGroup +
plan_self + plan_partner +
day +
# Random effects
(persuasion_self_cw + persuasion_partner_cw +
pressure_self_cw + pressure_partner_cw +
pushing_self_cw + pushing_partner_cw | coupleID)
, decomp = 'QR'
#, autocor = autocor_str
)
prior1 <- c(
brms::set_prior("normal(0, 5)", class = "b")
,brms::set_prior("normal(0, 20)", class = "Intercept", lb=1, ub=6)
, brms::set_prior("normal(0, 2)", class = "sd", group = "coupleID", lb = 0)
, brms::set_prior("student_t(3, 0, 2.5)", class = "sigma", lb = 0)
)
#brms::validate_prior(
# prior1,
# formula = formula,
# data = df_double,
# family = gaussian()
# )
#df_minimal <- df_double[, c("AorB", all.vars(as.formula(formula)))]
mood_gauss <- my_brm(
mi = use_mi,
imputed_data = implist,
formula = formula,
prior = prior1,
data = df_double,
family = gaussian(),
#control = list(adapt_delta = 0.95, max_treedepth = 15),
iter = iterations,
warmup = warmup,
chains = 4,
cores = 4,
seed = 7777,
file = file.path("models_cache_brms", paste0("mood_gauss", suffix))
)## Warning: Rows containing NAs were excluded from the model.
summary_mood <- summarize_brms(
mood_gauss,
stats_to_report = stats_to_report,
model_rows_fixed = model_rows_fixed,
model_rows_random = model_rows_random,
model_rownames_fixed = model_rownames_fixed,
model_rownames_random = model_rownames_random,
exponentiate = F) ## Warning in compute_rope(range = rope_range): Collinearity detected. Some
## VIFs are > 10. This may invalidate ROPE inferences!
Low Correlation
Term VIF VIF 95% CI Increased SE Tolerance
persuasion_self_cw 1.14 [ 1.11, 1.18] 1.07 0.88
persuasion_partner_cw 1.10 [ 1.07, 1.14] 1.05 0.91 pressure_self_cw 1.10 [ 1.07, 1.14] 1.05 0.91 pressure_partner_cw 1.07 [ 1.04, 1.11] 1.03 0.93 pushing_self_cw 1.16 [ 1.13, 1.21] 1.08 0.86 pushing_partner_cw 1.15 [ 1.11, 1.19] 1.07 0.87 isWeekend 1.04 [ 1.02, 1.09] 1.02 0.96 got_JITAI_self 1.09 [ 1.06, 1.13] 1.04 0.92 skilled_support 2.04 [ 1.96, 2.14] 1.43 0.49 studyGroup 1.13 [ 1.09, 1.17] 1.06 0.89 plan_self 1.29 [ 1.25, 1.34] 1.14 0.78 plan_partner 1.28 [ 1.23, 1.33] 1.13 0.78 day 1.99 [ 1.91, 2.08] 1.41 0.50 Tolerance 95% CI [0.85, 0.90] [0.88, 0.94] [0.88, 0.94] [0.90, 0.96] [0.83, 0.89] [0.84, 0.90] [0.92, 0.98] [0.89, 0.94] [0.47, 0.51] [0.86, 0.91] [0.75, 0.80] [0.75, 0.81] [0.48, 0.52]
Moderate Correlation
Term VIF VIF 95% CI Increased SE Tolerance
pressure_self_cb 5.52 [ 5.24, 5.83] 2.35 0.18
pressure_partner_cb 5.56 [ 5.27, 5.87] 2.36 0.18 pushing_self_cb 8.67 [ 8.20, 9.17] 2.94 0.12 pushing_partner_cb 8.62 [ 8.16, 9.12] 2.94 0.12 Tolerance 95% CI [0.17, 0.19] [0.17, 0.19] [0.11, 0.12] [0.11, 0.12]
High Correlation
Term VIF VIF 95% CI Increased SE Tolerance
persuasion_self_cb 13.89 [13.12, 14.71] 3.73 0.07
persuasion_partner_cb 13.87 [13.10, 14.68] 3.72 0.07 Tolerance 95% CI [0.07, 0.08] [0.07, 0.08]
## Sampling priors, please wait...
| Est. | SE | 95% CI | pd | ROPE | inside ROPE | BF | BF_Evidence | Rhat | Bulk_ESS | Tail_ESS | |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Intercept | 3.56*** | 0.16 | [ 3.23, 3.89] | 1.000 | [-0.11, 0.11] | 0.000 | >100 | Overwhelming Evidence | 1.000 | 5844 | 11807 |
| Within-Person Effects | |||||||||||
| Daily individual’s experienced persuasion | 0.00 | 0.02 | [-0.04, 0.04] | 0.522 | [-0.11, 0.11] | 1.000 | 0.004 | Very Strong Evidence for Null | 1.000 | 40415 | 29542 |
| Daily partner’s experienced persuasion | 0.02 | 0.02 | [-0.03, 0.06] | 0.757 | [-0.11, 0.11] | 1.000 | 0.005 | Very Strong Evidence for Null | 1.000 | 32705 | 29362 |
| Daily individual’s experienced pressure | -0.02 | 0.05 | [-0.14, 0.08] | 0.662 | [-0.11, 0.11] | 0.941 | 0.004 | Very Strong Evidence for Null | 1.000 | 39270 | 27760 |
| Daily partner’s experienced pressure | -0.01 | 0.05 | [-0.13, 0.10] | 0.584 | [-0.11, 0.11] | 0.945 | 0.004 | Very Strong Evidence for Null | 1.000 | 36823 | 27253 |
| Daily individual’s experienced pushing | 0.00 | 0.03 | [-0.06, 0.07] | 0.542 | [-0.11, 0.11] | 0.999 | 0.004 | Very Strong Evidence for Null | 1.000 | 45099 | 29788 |
| Daily partner’s experienced pushing | 0.07* | 0.03 | [ 0.00, 0.14] | 0.975 | [-0.11, 0.11] | 0.905 | 0.031 | Strong Evidence for Null | 1.000 | 33993 | 26876 |
| Day | 0.23** | 0.08 | [ 0.08, 0.38] | 0.998 | [-0.11, 0.11] | 0.072 | 0.246 | Moderate Evidence for Null | 1.000 | 59348 | 29338 |
| Own Actionplan | 0.13** | 0.04 | [ 0.05, 0.20] | 0.999 | [-0.11, 0.11] | 0.372 | 0.679 | Weak Evidence for Null | 1.000 | 61807 | 29764 |
| Partner Actionplan | -0.03 | 0.04 | [-0.10, 0.05] | 0.745 | [-0.11, 0.11] | 0.988 | 0.004 | Very Strong Evidence for Null | 1.000 | 58313 | 29688 |
| Daily wear time | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Daily support received | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Daily support provided (partner’s view) | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Is a Weekend | 0.15*** | 0.04 | [ 0.08, 0.21] | 1.000 | [-0.11, 0.11] | 0.193 | 4.224 | Moderate Evidence | 1.000 | 60906 | 27963 |
| JITAI received | -0.10* | 0.05 | [-0.20, -0.01] | 0.984 | [-0.11, 0.11] | 0.609 | 0.032 | Strong Evidence for Null | 1.000 | 81290 | 29772 |
| Days post skilled support intervention | 0.05 | 0.06 | [-0.06, 0.17] | 0.822 | [-0.11, 0.11] | 0.852 | 0.005 | Very Strong Evidence for Null | 1.000 | 47823 | 35802 |
| Between-Person Effects | |||||||||||
| Mean individual’s experienced persuasion | 0.39 | 0.27 | [-0.15, 0.94] | 0.925 | [-0.11, 0.11] | 0.122 | 0.036 | Strong Evidence for Null | 1.000 | 6491 | 11090 |
| Mean partner’s experienced persuasion | 0.26 | 0.27 | [-0.28, 0.81] | 0.833 | [-0.11, 0.11] | 0.208 | 0.021 | Very Strong Evidence for Null | 1.000 | 6363 | 11215 |
| Mean individual’s experienced pressure | -0.28 | 0.27 | [-0.81, 0.24] | 0.854 | [-0.11, 0.11] | 0.199 | 0.014 | Very Strong Evidence for Null | 1.000 | 7562 | 13793 |
| Mean partner’s experienced pressure | -0.29 | 0.26 | [-0.82, 0.24] | 0.863 | [-0.11, 0.11] | 0.189 | 0.015 | Very Strong Evidence for Null | 1.000 | 7293 | 14064 |
| Mean individual’s experienced pushing | 0.11 | 0.38 | [-0.66, 0.87] | 0.617 | [-0.11, 0.11] | 0.226 | 0.011 | Very Strong Evidence for Null | 1.001 | 8186 | 13586 |
| Mean partner’s experienced pushing | 0.27 | 0.38 | [-0.50, 1.02] | 0.764 | [-0.11, 0.11] | 0.185 | 0.013 | Very Strong Evidence for Null | 1.001 | 8168 | 13780 |
| Mean wear time | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Mean individual’s received support | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Mean partner’s received support | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Difference study group 2 | -0.21 | 0.22 | [-0.66, 0.24] | 0.830 | [-0.11, 0.11] | 0.258 | 0.028 | Very Strong Evidence for Null | 1.001 | 6328 | 10564 |
| Difference study group 3 | 0.35 | 0.22 | [-0.10, 0.80] | 0.937 | [-0.11, 0.11] | 0.126 | 0.055 | Strong Evidence for Null | 1.000 | 4708 | 9432 |
| Random Effects | |||||||||||
| sd(Intercept) | 0.56 | 0.07 | [0.44, 0.74] | NA | NA | NA | NA | NA | 1.000 | 9471 | 16569 |
| sd(Daily individual’s experienced persuasion) | 0.04 | 0.03 | [0.00, 0.10] | NA | NA | NA | NA | NA | 1.000 | 11769 | 14579 |
| sd(Daily partner’s experienced persuasion) | 0.07 | 0.03 | [0.01, 0.13] | NA | NA | NA | NA | NA | 1.000 | 9695 | 8040 |
| sd(Daily individual’s experienced pressure) | 0.07 | 0.06 | [0.00, 0.25] | NA | NA | NA | NA | NA | 1.000 | 16105 | 17953 |
| sd(Daily partner’s experienced pressure) | 0.08 | 0.07 | [0.00, 0.27] | NA | NA | NA | NA | NA | 1.000 | 14513 | 16913 |
| sd(Daily individual’s experienced pushing) | 0.05 | 0.04 | [0.00, 0.14] | NA | NA | NA | NA | NA | 1.000 | 15105 | 17267 |
| sd(Daily partner’s experienced pushing) | 0.06 | 0.05 | [0.00, 0.16] | NA | NA | NA | NA | NA | 1.000 | 13248 | 15418 |
| Additional Parameters | |||||||||||
| sigma | 0.96 | 0.01 | [0.94, 0.98] | NA | NA | NA | NA | NA | 1.000 | 59952 | 29269 |
plot(
bayestestR::p_direction(mood_gauss),
priors = TRUE
) +
coord_cartesian(xlim = c(-3, 3)) +
theme_bw()## Warning in `==.default`(dens$Parameter, parameter): longer object length
## is not a multiple of shorter object length
## Warning in is.na(e1) | is.na(e2): longer object length is not a multiple
## of shorter object length
## Possible multicollinearity between b_pressure_self_cb and
## b_persuasion_self_cb (r = 0.8), b_pressure_partner_cb and
## b_persuasion_self_cb (r = 0.77), b_pressure_self_cb and
## b_persuasion_partner_cb (r = 0.77), b_pressure_partner_cb and
## b_persuasion_partner_cb (r = 0.8), b_pressure_partner_cb and
## b_pressure_self_cb (r = 0.75), b_pushing_partner_cb and
## b_pushing_self_cb (r = 0.89). This might lead to inappropriate
## results. See 'Details' in '?rope'.
conditional_spaghetti(
mood_gauss,
effects = c('pushing_partner_cw'),
group_var = 'coupleID',
plot_full_range = TRUE
)$pushing_partner_cw
df_double$reactance_ordinal <- factor(df_double$reactance,
levels = 0:5,
ordered = TRUE)
formula <- bf(
reactance_ordinal ~
persuasion_self_cw + persuasion_partner_cw +
pressure_self_cw + pressure_partner_cw +
pushing_self_cw + pushing_partner_cw +
#support_self_cw + support_partner_cw +
isWeekend +
got_JITAI_self +
skilled_support +
persuasion_self_cb + persuasion_partner_cb +
pressure_self_cb + pressure_partner_cb +
pushing_self_cb + pushing_partner_cb +
#support_self_cb + support_partner_cb +
studyGroup +
plan_self + plan_partner +
day +
# Random effects
(persuasion_self_cw + persuasion_partner_cw +
pressure_self_cw + pressure_partner_cw +
pushing_self_cw + pushing_partner_cw | coupleID)
, decomp = 'QR'
#, autocor = autocor_str
)
prior1 <- c(
brms::set_prior("normal(0, 2.5)", class = "b")
, brms::set_prior("normal(0, 2)", class = "sd", group = "coupleID", lb = 0)
)
#brms::validate_prior(
# prior1,
# formula = formula,
# data = df_double,
# family = cumulative() # HURDLE_CUMULATIVE
# )
#df_minimal <- df_double[, c("AorB", all.vars(as.formula(formula)))]
reactance_ordinal <- my_brm(
mi = use_mi,
imputed_data = implist,
formula = formula,
prior = prior1,
data = df_double,
family = brms::cumulative(),
#control = list(adapt_delta = 0.95),
iter = iterations,
warmup = warmup,
chains = 4,
cores = 4,
seed = 7777
, file = file.path("models_cache_brms", paste0("reactance_ordinal", suffix))
)## Warning: Rows containing NAs were excluded from the model.
summary_reactance_ordinal <- summarize_brms(
reactance_ordinal,
stats_to_report = stats_to_report,
rope_range = c(-0.18, 0.18),
model_rows_fixed = model_rows_fixed_ordinal,
model_rows_random = model_rows_random_ordinal,
model_rownames_fixed = model_rownames_fixed_ordinal,
model_rownames_random = model_rownames_random_ordinal,
exponentiate = T) ## Warning in compute_rope(range = rope_range): Collinearity detected. Some
## VIFs are > 10. This may invalidate ROPE inferences!
Low Correlation
Term VIF VIF 95% CI Increased SE Tolerance
pressure_partner_cw 1.72 [ 1.63, 1.83] 1.31 0.58 pushing_self_cw 1.33 [ 1.27, 1.41] 1.15 0.75 pushing_partner_cw 1.11 [ 1.07, 1.18] 1.06 0.90 isWeekend 1.10 [ 1.06, 1.16] 1.05 0.91 got_JITAI_self 1.08 [ 1.04, 1.14] 1.04 0.93 skilled_support 1.39 [ 1.33, 1.48] 1.18 0.72 persuasion_self_cb 1.20 [ 1.15, 1.26] 1.09 0.83 persuasion_partner_cb 1.15 [ 1.11, 1.22] 1.07 0.87 pressure_self_cb 1.26 [ 1.20, 1.33] 1.12 0.79 pressure_partner_cb 4.11 [ 3.83, 4.42] 2.03 0.24 plan_self 2.34 [ 2.20, 2.50] 1.53 0.43 plan_partner 2.66 [ 2.49, 2.85] 1.63 0.38 day 2.23 [ 2.09, 2.38] 1.49 0.45 Tolerance 95% CI [0.55, 0.61] [0.71, 0.79] [0.85, 0.93] [0.86, 0.95] [0.88, 0.96] [0.68, 0.75] [0.79, 0.87] [0.82, 0.90] [0.75, 0.83] [0.23, 0.26] [0.40, 0.46] [0.35, 0.40] [0.42, 0.48]
Moderate Correlation
Term VIF VIF 95% CI Increased SE Tolerance
pressure_self_cw 6.37 [ 5.92, 6.87] 2.52 0.16 pushing_self_cb 6.52 [ 6.05, 7.04] 2.55 0.15 pushing_partner_cb 8.27 [ 7.66, 8.93] 2.88 0.12 Tolerance 95% CI [0.15, 0.17] [0.14, 0.17] [0.11, 0.13]
High Correlation
Term VIF VIF 95% CI Increased SE Tolerance
persuasion_self_cw 14.96 [13.83, 16.19] 3.87 0.07
persuasion_partner_cw 15.75 [14.56, 17.05] 3.97 0.06 studyGroup 22.18 [20.49, 24.03] 4.71 0.05 Tolerance 95% CI [0.06, 0.07] [0.06, 0.07] [0.04, 0.05]
## Sampling priors, please wait...
===> warning: knot 9 removed - double knot ===> warning: knot 10 removed - double knot ===> warning: knot 11 removed - double knot ===> warning: knot 12 removed - double knot ===> warning: knot 13 removed - double knot ===> warning: knot 14 removed - double knot ===> warning: knot 15 removed - double knot ===> warning: knot 16 removed - double knot ===> warning: knot 3 removed - double knot ===> warning: knot 4 removed - double knot ===> warning: knot 5 removed - double knot ===> warning: knot 6 removed - double knot ===> warning: knot 7 removed - double knot ===> warning: knot 8 removed - double knot * several double knots suggests that your data is strongly rounded, attention might be required ===> warning: knot 9 removed - double knot ===> warning: knot 10 removed - double knot ===> warning: knot 11 removed - double knot ===> warning: knot 12 removed - double knot ===> warning: knot 13 removed - double knot ===> warning: knot 14 removed - double knot ===> warning: knot 15 removed - double knot ===> warning: knot 16 removed - double knot ===> warning: knot 3 removed - double knot ===> warning: knot 4 removed - double knot ===> warning: knot 5 removed - double knot ===> warning: knot 6 removed - double knot ===> warning: knot 7 removed - double knot ===> warning: knot 8 removed - double knot several double knots suggests that your data is strongly rounded, attention might be required *
| OR | SE | 95% CI | pd | ROPE | inside ROPE | BF | BF_Evidence | Rhat | Bulk_ESS | Tail_ESS | |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Intercepts | |||||||||||
| Intercept | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Intercept[1] | 2.00 | 0.79 | [ 0.93, 4.49] | 0.963 | [0.84, 1.20] | 0.080 | 0.176 | Moderate Evidence for Null | 1.000 | 30668 | 29360 |
| Intercept[2] | 4.41*** | 1.79 | [ 2.05, 10.05] | 1.000 | [0.84, 1.20] | 0.000 | 26.556 | Strong Evidence | 1.000 | 30524 | 27684 |
| Intercept[3] | 11.92*** | 4.98 | [ 5.40, 28.22] | 1.000 | [0.84, 1.20] | 0.000 | >100 | Overwhelming Evidence | 1.000 | 30791 | 29212 |
| Intercept[4] | 50.98*** | 23.49 | [ 21.54, 131.03] | 1.000 | [0.84, 1.20] | 0.000 | >100 | Overwhelming Evidence | 1.000 | 30657 | 29001 |
| Intercept[5] | 1584.75*** | 1122.79 | [435.24, 6903.06] | 1.000 | [0.84, 1.20] | 0.000 | >100 | Overwhelming Evidence | 1.000 | 36368 | 30707 |
| Within-Person Effects | |||||||||||
| Daily individual’s experienced persuasion | 0.86 | 0.07 | [ 0.72, 1.01] | 0.966 | [0.84, 1.20] | 0.613 | 0.248 | Moderate Evidence for Null | 1.000 | 30083 | 26859 |
| Daily partner’s experienced persuasion | 1.01 | 0.10 | [ 0.82, 1.24] | 0.554 | [0.84, 1.20] | 0.916 | 0.041 | Strong Evidence for Null | 1.000 | 29423 | 25629 |
| Daily individual’s experienced pressure | 1.83* | 0.35 | [ 1.15, 2.62] | 0.993 | [0.84, 1.20] | 0.032 | 1.803 | Weak Evidence | 1.000 | 17725 | 19483 |
| Daily partner’s experienced pressure | 1.23 | 0.30 | [ 0.68, 2.10] | 0.797 | [0.84, 1.20] | 0.382 | 0.083 | Strong Evidence for Null | 1.000 | 18931 | 16192 |
| Daily individual’s experienced pushing | 1.24 | 0.14 | [ 0.99, 1.56] | 0.967 | [0.84, 1.20] | 0.390 | 0.229 | Moderate Evidence for Null | 1.000 | 26253 | 25762 |
| Daily partner’s experienced pushing | 0.96 | 0.13 | [ 0.72, 1.25] | 0.625 | [0.84, 1.20] | 0.795 | 0.044 | Strong Evidence for Null | 1.000 | 26476 | 23939 |
| Day | 1.72 | 0.82 | [ 0.66, 4.37] | 0.870 | [0.84, 1.20] | 0.158 | 0.066 | Strong Evidence for Null | 1.000 | 38344 | 29944 |
| Own Actionplan | 0.79 | 0.24 | [ 0.44, 1.41] | 0.782 | [0.84, 1.20] | 0.351 | 0.052 | Strong Evidence for Null | 1.000 | 39810 | 29550 |
| Partner Actionplan | 0.86 | 0.23 | [ 0.52, 1.45] | 0.712 | [0.84, 1.20] | 0.445 | 0.043 | Strong Evidence for Null | 1.000 | 36678 | 30639 |
| Daily wear time | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Daily support received | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Daily support provided (partner’s view) | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Is a Weekend | 0.62* | 0.13 | [ 0.40, 0.94] | 0.988 | [0.84, 1.20] | 0.073 | 0.435 | Weak Evidence for Null | 1.000 | 47237 | 32149 |
| JITAI received | 1.27 | 0.35 | [ 0.74, 2.18] | 0.813 | [0.84, 1.20] | 0.350 | 0.054 | Strong Evidence for Null | 1.000 | 46810 | 31201 |
| Days post skilled support intervention | 0.82 | 0.29 | [ 0.41, 1.65] | 0.712 | [0.84, 1.20] | 0.336 | 0.042 | Strong Evidence for Null | 1.000 | 47691 | 29126 |
| Between-Person Effects | |||||||||||
| Mean individual’s experienced persuasion | 0.87 | 0.43 | [ 0.32, 2.38] | 0.613 | [0.84, 1.20] | 0.273 | 0.054 | Strong Evidence for Null | 1.000 | 21453 | 25466 |
| Mean partner’s experienced persuasion | 1.33 | 0.77 | [ 0.44, 4.30] | 0.695 | [0.84, 1.20] | 0.219 | 0.060 | Strong Evidence for Null | 1.000 | 21944 | 24421 |
| Mean individual’s experienced pressure | 3.91* | 2.12 | [ 1.38, 12.13] | 0.994 | [0.84, 1.20] | 0.011 | 1.192 | Weak Evidence | 1.000 | 19938 | 25916 |
| Mean partner’s experienced pressure | 1.29 | 0.74 | [ 0.40, 3.90] | 0.675 | [0.84, 1.20] | 0.223 | 0.052 | Strong Evidence for Null | 1.000 | 18928 | 24521 |
| Mean individual’s experienced pushing | 1.54 | 1.14 | [ 0.36, 6.99] | 0.723 | [0.84, 1.20] | 0.165 | 0.061 | Strong Evidence for Null | 1.000 | 17717 | 22675 |
| Mean partner’s experienced pushing | 0.08** | 0.08 | [ 0.01, 0.52] | 0.996 | [0.84, 1.20] | 0.005 | 2.076 | Weak Evidence | 1.000 | 21162 | 25207 |
| Mean wear time | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Mean individual’s received support | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Mean partner’s received support | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Difference study group 2 | 1.08 | 0.41 | [ 0.49, 2.32] | 0.578 | [0.84, 1.20] | 0.356 | 0.056 | Strong Evidence for Null | 1.000 | 20601 | 26250 |
| Difference study group 3 | 0.36* | 0.15 | [ 0.15, 0.87] | 0.987 | [0.84, 1.20] | 0.025 | 0.836 | Weak Evidence for Null | 1.000 | 18739 | 24668 |
| Random Effects | |||||||||||
| sd(Intercept) | 0.64 | 0.22 | [0.21, 1.14] | NA | NA | NA | NA | NA | 1.000 | 8500 | 8003 |
| sd(Daily individual’s experienced persuasion) | 0.18 | 0.13 | [0.01, 0.45] | NA | NA | NA | NA | NA | 1.000 | 6069 | 10323 |
| sd(Daily partner’s experienced persuasion) | 0.22 | 0.14 | [0.01, 0.53] | NA | NA | NA | NA | NA | 1.000 | 9304 | 11594 |
| sd(Daily individual’s experienced pressure) | 0.53 | 0.27 | [0.05, 1.14] | NA | NA | NA | NA | NA | 1.001 | 7911 | 9419 |
| sd(Daily partner’s experienced pressure) | 0.41 | 0.38 | [0.02, 1.55] | NA | NA | NA | NA | NA | 1.001 | 9709 | 15866 |
| sd(Daily individual’s experienced pushing) | 0.24 | 0.14 | [0.02, 0.56] | NA | NA | NA | NA | NA | 1.000 | 8773 | 11200 |
| sd(Daily partner’s experienced pushing) | 0.15 | 0.14 | [0.01, 0.61] | NA | NA | NA | NA | NA | 1.000 | 13960 | 16216 |
| Additional Parameters | |||||||||||
| sigma | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| disc | 1.00 | 0.00 | [1.00, 1.00] | NA | NA | NA | NA | NA | NA | NA | NA |
plot(
bayestestR::p_direction(reactance_ordinal),
priors = TRUE
) +
coord_cartesian(xlim = c(-6, 6)) +
theme_bw()## Possible multicollinearity between b_Intercept[4] and b_Intercept[2]
## (r = 0.87), b_pressure_partner_cb and b_persuasion_partner_cb (r =
## 0.78). This might lead to inappropriate results. See 'Details' in
## '?rope'.
introduce_binary_reactance <- function(data) {
data$is_reactance <- factor(data$reactance > 0, levels = c(FALSE, TRUE), labels = c(0, 1))
return(data)
}
df_double <- introduce_binary_reactance(df_double)
if (use_mi) {
for (i in seq_along(implist)) {
implist[[i]] <- introduce_binary_reactance(implist[[i]])
}
}
formula <- bf(
is_reactance ~
persuasion_self_cw + persuasion_partner_cw +
pressure_self_cw + pressure_partner_cw +
pushing_self_cw + pushing_partner_cw +
#support_self_cw + support_partner_cw +
isWeekend +
got_JITAI_self +
skilled_support +
persuasion_self_cb + persuasion_partner_cb +
pressure_self_cb + pressure_partner_cb +
pushing_self_cb + pushing_partner_cb +
#support_self_cb + support_partner_cb +
studyGroup +
plan_self + plan_partner +
day +
# Random effects
(persuasion_self_cw + persuasion_partner_cw +
pressure_self_cw + pressure_partner_cw +
pushing_self_cw + pushing_partner_cw | coupleID)
, decomp = 'QR'
#, autocor = autocor_str
)
prior1 <- c(
brms::set_prior("normal(0, 2.5)", class = "b")
, brms::set_prior("normal(0, 10)", class = "Intercept", lb=0, ub=5)
, brms::set_prior("normal(0, 2)", class = "sd", group = "coupleID", lb = 0)
)
#brms::validate_prior(
# prior1,
# formula = formula,
# data = df_double,
# family = bernoulli()
# )
#df_minimal <- df_double[, c("AorB", all.vars(as.formula(formula)))]
is_reactance <- my_brm(
mi = use_mi,
imputed_data = implist,
formula = formula,
prior = prior1,
data = df_double,
family = brms::bernoulli(),
#control = list(adapt_delta = 0.95, max_treedepth = 15),
iter = iterations,
warmup = warmup,
chains = 4,
cores = 4,
seed = 7777,
file = file.path("models_cache_brms", paste0("is_reactance", suffix))
#, file_refit = 'always'
)## Warning: Rows containing NAs were excluded from the model.
summary_is_reactance <- summarize_brms(
is_reactance,
stats_to_report = stats_to_report,
model_rows_fixed = model_rows_fixed,
model_rows_random = model_rows_random,
model_rownames_fixed = model_rownames_fixed,
model_rownames_random = model_rownames_random,
exponentiate = T) ## Sampling priors, please wait...
| OR | SE | 95% CI | pd | ROPE | inside ROPE | BF | BF_Evidence | Rhat | Bulk_ESS | Tail_ESS | |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Intercept | 0.62 | 0.30 | [0.23, 1.62] | 0.837 | [0.83, 1.20] | 0.184 | 0.085 | Strong Evidence for Null | 1.000 | 25105 | 29851 |
| Within-Person Effects | |||||||||||
| Daily individual’s experienced persuasion | 0.84 | 0.09 | [0.68, 1.03] | 0.958 | [0.83, 1.20] | 0.541 | 0.242 | Moderate Evidence for Null | 1.000 | 31253 | 28630 |
| Daily partner’s experienced persuasion | 1.14 | 0.18 | [0.84, 1.60] | 0.800 | [0.83, 1.20] | 0.612 | 0.095 | Strong Evidence for Null | 1.000 | 19207 | 24180 |
| Daily individual’s experienced pressure | 1.99* | 0.66 | [1.03, 4.87] | 0.979 | [0.83, 1.20] | 0.051 | 0.958 | Weak Evidence for Null | 1.001 | 16301 | 17024 |
| Daily partner’s experienced pressure | 1.45 | 0.65 | [0.56, 4.71] | 0.802 | [0.83, 1.20] | 0.223 | 0.137 | Moderate Evidence for Null | 1.000 | 17230 | 18755 |
| Daily individual’s experienced pushing | 1.37* | 0.19 | [1.06, 1.83] | 0.991 | [0.83, 1.20] | 0.153 | 1.004 | Weak Evidence | 1.000 | 29562 | 27215 |
| Daily partner’s experienced pushing | 0.91 | 0.18 | [0.60, 1.36] | 0.687 | [0.83, 1.20] | 0.577 | 0.080 | Strong Evidence for Null | 1.000 | 26807 | 25067 |
| Day | 2.11 | 1.20 | [0.69, 6.36] | 0.906 | [0.83, 1.20] | 0.108 | 0.105 | Moderate Evidence for Null | 1.000 | 38208 | 29162 |
| Own Actionplan | 0.81 | 0.27 | [0.42, 1.57] | 0.737 | [0.83, 1.20] | 0.338 | 0.055 | Strong Evidence for Null | 1.000 | 40855 | 28916 |
| Partner Actionplan | 0.84 | 0.25 | [0.46, 1.52] | 0.719 | [0.83, 1.20] | 0.386 | 0.052 | Strong Evidence for Null | 1.000 | 40714 | 27635 |
| Daily wear time | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Daily support received | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Daily support provided (partner’s view) | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Is a Weekend | 0.61 | 0.16 | [0.36, 1.02] | 0.970 | [0.83, 1.20] | 0.116 | 0.260 | Moderate Evidence for Null | 1.000 | 48766 | 31499 |
| JITAI received | 1.47 | 0.48 | [0.77, 2.78] | 0.877 | [0.83, 1.20] | 0.225 | 0.085 | Strong Evidence for Null | 1.000 | 40721 | 30512 |
| Days post skilled support intervention | 0.72 | 0.29 | [0.33, 1.56] | 0.795 | [0.83, 1.20] | 0.254 | 0.062 | Strong Evidence for Null | 1.000 | 42039 | 33004 |
| Between-Person Effects | |||||||||||
| Mean individual’s experienced persuasion | 1.59 | 1.09 | [0.42, 6.85] | 0.751 | [0.83, 1.20] | 0.170 | 0.089 | Strong Evidence for Null | 1.000 | 18671 | 26193 |
| Mean partner’s experienced persuasion | 1.97 | 1.49 | [0.46, 9.66] | 0.820 | [0.83, 1.20] | 0.132 | 0.103 | Moderate Evidence for Null | 1.000 | 18982 | 26195 |
| Mean individual’s experienced pressure | 36.37** | 48.21 | [3.22, 617.66] | 0.999 | [0.83, 1.20] | 0.002 | 9.372 | Moderate Evidence | 1.000 | 12515 | 18763 |
| Mean partner’s experienced pressure | 2.00 | 2.78 | [0.10, 24.54] | 0.689 | [0.83, 1.20] | 0.089 | 0.136 | Moderate Evidence for Null | 1.000 | 11967 | 21167 |
| Mean individual’s experienced pushing | 0.73 | 0.86 | [0.07, 7.13] | 0.604 | [0.83, 1.20] | 0.119 | 0.084 | Strong Evidence for Null | 1.000 | 16764 | 23152 |
| Mean partner’s experienced pushing | 0.05* | 0.07 | [0.00, 0.64] | 0.989 | [0.83, 1.20] | 0.008 | 1.171 | Weak Evidence | 1.000 | 16820 | 23993 |
| Mean wear time | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Mean individual’s received support | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Mean partner’s received support | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Difference study group 2 | 1.10 | 0.60 | [0.36, 3.27] | 0.570 | [0.83, 1.20] | 0.257 | 0.080 | Strong Evidence for Null | 1.000 | 15870 | 22252 |
| Difference study group 3 | 0.36 | 0.21 | [0.11, 1.21] | 0.955 | [0.83, 1.20] | 0.053 | 0.363 | Weak Evidence for Null | 1.000 | 16863 | 22561 |
| Random Effects | |||||||||||
| sd(Intercept) | 1.17 | 0.28 | [0.69, 1.83] | NA | NA | NA | NA | NA | 1.001 | 11762 | 18857 |
| sd(Daily individual’s experienced persuasion) | 0.25 | 0.15 | [0.02, 0.57] | NA | NA | NA | NA | NA | 1.001 | 6791 | 11619 |
| sd(Daily partner’s experienced persuasion) | 0.52 | 0.21 | [0.14, 1.03] | NA | NA | NA | NA | NA | 1.000 | 10860 | 10273 |
| sd(Daily individual’s experienced pressure) | 1.10 | 0.59 | [0.11, 2.49] | NA | NA | NA | NA | NA | 1.001 | 6782 | 6614 |
| sd(Daily partner’s experienced pressure) | 0.87 | 0.71 | [0.04, 2.83] | NA | NA | NA | NA | NA | 1.000 | 11338 | 14360 |
| sd(Daily individual’s experienced pushing) | 0.26 | 0.17 | [0.02, 0.65] | NA | NA | NA | NA | NA | 1.000 | 10950 | 13641 |
| sd(Daily partner’s experienced pushing) | 0.25 | 0.22 | [0.01, 0.94] | NA | NA | NA | NA | NA | 1.000 | 14970 | 18666 |
| Additional Parameters | |||||||||||
| sigma | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
plot(
bayestestR::p_direction(is_reactance),
priors = TRUE
) +
coord_cartesian(xlim = c(-6, 6)) +
theme_bw()## Warning in `==.default`(dens$Parameter, parameter): longer object length
## is not a multiple of shorter object length
## Warning in is.na(e1) | is.na(e2): longer object length is not a multiple
## of shorter object length
conditional_spaghetti(
is_reactance,
effects = c('pressure_self_cw', 'pushing_self_cw'),
group_var = 'coupleID',
plot_full_range = TRUE
)\(pressure_self_cw
<img
src="04_SensitivityOtherCovariates_files/figure-html/report_is_reactance-3.png"
width="2400" />\)pushing_self_cw
summary_all_models <- report_side_by_side(
pa_sub,
pa_obj_log,
mood_gauss,
is_reactance,
stats_to_report = c('CI'),
model_rows_random = model_rows_random,
model_rows_fixed = model_rows_fixed,
model_rownames_random = model_rownames_random,
model_rownames_fixed = model_rownames_fixed
)[1] “pa_sub”
## Warning in summarize_brms(model, exponentiate = exponentiate,
## stats_to_report = stats_to_report, : Coefficients were exponentiated.
## Double check if this was intended.
[1] “pa_obj_log”
## Warning in summarize_brms(model, exponentiate = exponentiate,
## stats_to_report = stats_to_report, : Coefficients were exponentiated.
## Double check if this was intended.
[1] “mood_gauss” [1] “is_reactance”
summary_all_models <- summary_all_models %>%
print_df(rows_to_pack = rows_to_pack) %>%
add_header_above(
c(
" ", "Hurdle Component" = 2, "Non-Zero Component" = 2,
" " = 6
)
) %>%
add_header_above(
c(" ", "Subjective MVPA Hurdle Lognormal" = 4,
"Device-Based MVPA Log (Gaussian)" = 2,
"Mood Gaussian" = 2,
#"Reactance Ordinal" = 2,
"Reactance Dichotome" = 2
)
)
export_xlsx(
summary_all_models,
rows_to_pack = rows_to_pack,
file.path("Output", paste0("AllModels", suffix, ".xlsx")),
merge_option = 'both',
simplify_2nd_row = TRUE,
line_above_rows = c(1,2),
line_below_rows = c(-1)
)
summary_all_models| exp(Est.)_hu pa_sub | 95% CI_hu pa_sub | exp(Est.)_nonzero pa_sub | 95% CI_nonzero pa_sub | exp(Est.) pa_obj_log | 95% CI pa_obj_log | Est. mood_gauss | 95% CI mood_gauss | OR is_reactance | 95% CI is_reactance | |
|---|---|---|---|---|---|---|---|---|---|---|
| Intercept | 0.18*** | [ 0.10, 0.30] | 29.02*** | [23.19, 36.50] | 108.21*** | [90.06, 129.44] | 3.56*** | [ 3.23, 3.89] | 0.62 | [0.23, 1.62] |
| Within-Person Effects | ||||||||||
| Daily individual’s experienced persuasion | 1.57*** | [ 1.38, 1.83] | 1.03 | [ 0.98, 1.08] | 1.03 | [ 1.00, 1.06] | 0.00 | [-0.04, 0.04] | 0.84 | [0.68, 1.03] |
| Daily partner’s experienced persuasion | 1.33*** | [ 1.18, 1.53] | 1.03 | [ 0.98, 1.07] | 1.01 | [ 0.98, 1.05] | 0.02 | [-0.03, 0.06] | 1.14 | [0.84, 1.60] |
| Daily individual’s experienced pressure | 0.98 | [ 0.70, 1.37] | 0.90* | [ 0.81, 1.00] | 0.95 | [ 0.88, 1.01] | -0.02 | [-0.14, 0.08] | 1.99* | [1.03, 4.87] |
| Daily partner’s experienced pressure | 1.56* | [ 1.04, 2.74] | 0.95 | [ 0.86, 1.04] | 0.98 | [ 0.92, 1.05] | -0.01 | [-0.13, 0.10] | 1.45 | [0.56, 4.71] |
| Daily individual’s experienced pushing | 1.00 | [ 0.76, 1.35] | 0.99 | [ 0.92, 1.05] | 1.02 | [ 0.96, 1.07] | 0.00 | [-0.06, 0.07] | 1.37* | [1.06, 1.83] |
| Daily partner’s experienced pushing | 1.33* | [ 1.05, 1.72] | 0.96 | [ 0.90, 1.02] | 1.00 | [ 0.96, 1.04] | 0.07* | [ 0.00, 0.14] | 0.91 | [0.60, 1.36] |
| Day | 0.69 | [ 0.46, 1.04] | 0.94 | [ 0.79, 1.11] | 0.94 | [ 0.86, 1.03] | 0.23** | [ 0.08, 0.38] | 2.11 | [0.69, 6.36] |
| Own Actionplan | 11.49*** | [ 9.35, 14.19] | 1.37*** | [ 1.24, 1.50] | 1.08** | [ 1.02, 1.13] | 0.13** | [ 0.05, 0.20] | 0.81 | [0.42, 1.57] |
| Partner Actionplan | 1.24* | [ 1.01, 1.51] | 1.08 | [ 0.99, 1.18] | 1.05* | [ 1.00, 1.11] | -0.03 | [-0.10, 0.05] | 0.84 | [0.46, 1.52] |
| Daily wear time | NA | NA | NA | NA | 1.00*** | [ 1.00, 1.00] | NA | NA | NA | NA |
| Daily support received | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Daily support provided (partner’s view) | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Is a Weekend | 1.57*** | [ 1.29, 1.91] | 1.19*** | [ 1.10, 1.28] | 1.06** | [ 1.01, 1.11] | 0.15*** | [ 0.08, 0.21] | 0.61 | [0.36, 1.02] |
| JITAI received | 0.49*** | [ 0.38, 0.63] | 0.94 | [ 0.85, 1.04] | 0.93* | [ 0.87, 0.99] | -0.10* | [-0.20, -0.01] | 1.47 | [0.77, 2.78] |
| Days post skilled support intervention | 1.31 | [ 0.96, 1.80] | 1.07 | [ 0.94, 1.21] | 1.05 | [ 0.98, 1.13] | 0.05 | [-0.06, 0.17] | 0.72 | [0.33, 1.56] |
| Between-Person Effects | ||||||||||
| Mean individual’s experienced persuasion | 1.44 | [ 0.65, 3.13] | 0.99 | [ 0.72, 1.34] | 1.10 | [ 0.82, 1.49] | 0.39 | [-0.15, 0.94] | 1.59 | [0.42, 6.85] |
| Mean partner’s experienced persuasion | 1.32 | [ 0.60, 2.88] | 0.93 | [ 0.68, 1.26] | 0.98 | [ 0.73, 1.33] | 0.26 | [-0.28, 0.81] | 1.97 | [0.46, 9.66] |
| Mean individual’s experienced pressure | 0.16*** | [ 0.05, 0.47] | 0.69 | [ 0.33, 1.45] | 0.98 | [ 0.72, 1.32] | -0.28 | [-0.81, 0.24] | 36.37** | [3.22, 617.66] |
| Mean partner’s experienced pressure | 0.27* | [ 0.08, 0.78] | 0.53 | [ 0.25, 1.12] | 0.98 | [ 0.73, 1.32] | -0.29 | [-0.82, 0.24] | 2.00 | [0.10, 24.54] |
| Mean individual’s experienced pushing | 1.39 | [ 0.42, 4.73] | 1.47 | [ 0.86, 2.52] | 0.96 | [ 0.63, 1.47] | 0.11 | [-0.66, 0.87] | 0.73 | [0.07, 7.13] |
| Mean partner’s experienced pushing | 2.33 | [ 0.70, 7.85] | 1.63 | [ 0.94, 2.83] | 1.20 | [ 0.79, 1.84] | 0.27 | [-0.50, 1.02] | 0.05* | [0.00, 0.64] |
| Mean wear time | NA | NA | NA | NA | 1.00 | [ 1.00, 1.00] | NA | NA | NA | NA |
| Mean individual’s received support | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Mean partner’s received support | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| Difference study group 2 | 0.79 | [ 0.39, 1.60] | 1.06 | [ 0.81, 1.38] | 0.94 | [ 0.73, 1.20] | -0.21 | [-0.66, 0.24] | 1.10 | [0.36, 3.27] |
| Difference study group 3 | 1.17 | [ 0.59, 2.33] | 1.25 | [ 0.97, 1.60] | 1.03 | [ 0.80, 1.33] | 0.35 | [-0.10, 0.80] | 0.36 | [0.11, 1.21] |
| Random Effects | ||||||||||
| sd(Intercept) | 0.83 | [0.63, 1.12] | 0.29 | [0.22, 0.40] | 0.30 | [0.24, 0.40] | 0.56 | [0.44, 0.74] | 1.17 | [0.69, 1.83] |
| sd(Daily individual’s experienced persuasion) | 0.20 | [0.03, 0.40] | 0.12 | [0.08, 0.17] | 0.05 | [0.02, 0.08] | 0.04 | [0.00, 0.10] | 0.25 | [0.02, 0.57] |
| sd(Daily partner’s experienced persuasion) | 0.16 | [0.01, 0.36] | 0.08 | [0.04, 0.13] | 0.05 | [0.02, 0.09] | 0.07 | [0.01, 0.13] | 0.52 | [0.14, 1.03] |
| sd(Daily individual’s experienced pressure) | 0.21 | [0.01, 0.78] | 0.07 | [0.00, 0.23] | 0.04 | [0.00, 0.13] | 0.07 | [0.00, 0.25] | 1.10 | [0.11, 2.49] |
| sd(Daily partner’s experienced pressure) | 0.39 | [0.02, 1.25] | 0.05 | [0.00, 0.18] | 0.03 | [0.00, 0.11] | 0.08 | [0.00, 0.27] | 0.87 | [0.04, 2.83] |
| sd(Daily individual’s experienced pushing) | 0.51 | [0.21, 0.91] | 0.08 | [0.01, 0.16] | 0.08 | [0.01, 0.16] | 0.05 | [0.00, 0.14] | 0.26 | [0.02, 0.65] |
| sd(Daily partner’s experienced pushing) | 0.23 | [0.01, 0.59] | 0.07 | [0.01, 0.14] | 0.03 | [0.00, 0.09] | 0.06 | [0.00, 0.16] | 0.25 | [0.01, 0.94] |
| Additional Parameters | ||||||||||
| sigma | NA | NA | 0.67 | [0.65, 0.70] | 0.57 | [0.56, 0.59] | 0.96 | [0.94, 0.98] | NA | NA |
Analyses were conducted using the R Statistical language (version 4.4.2; R Core Team, 2024) on Windows 11 x64 (build 26100)